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The invention concerns an image processing apparatus and a method for 
the reconstmction of time-dependent representations of an object, particularly during 
perfusion imaging. Moreover the invention relates to an X-ray examination system 
comprising such an apparatus and a computer program implementing said method as 

1 0 well as to an X-ray system. 

Perfusion imaging is an important functional imaging method in medical 
applications which delivers information about the blood supply of the tissue and the 
tissue viability as well as valuable information during interventions. According to 
cuirent techniques, perfusion imaging requires the acquisition of a time series (with 

1 5 typically about 40 members) of (multiple) cross-sectional images with a frame rate in 
the order of 0.75 - 1 second. After acquisition and reconstruction of the mdividual 
cross-sectional images, the image time series is evaluated and a perfusion map is 
computed. 

In a perfusion imaging procedure described in the US 6 373 920 Bl, a 
20 first set of projections is taken without contrast agent and later a second set of 
projections is taken with contrast agent. From the difference of the images 
reconstructed based on each set, a perfusion image is calculated. 

As current C-arm X-ray systems are limited in their rotational speed and 
detector frame rate, such systems cannot be used for perfusion imaging based on 
25 cross-sectional images generated in short time intervals. 

Based on this situation it was an object of the present invention to 
provide means for a facilitated perfusion imagmg, especially by use of normal C-arm 
based X-ray systems. 

30 This object is achieved by an image processing apparatus according to 

claim 1 , an X-ray examination system according to claim 9, a method according to 
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claim 13, a computer program according to claim 15, a record carrier according to 
claim 16, and an X-ray system according to Claim 17, Preferred embodiments are 

subject of the dependent claims. 

According to a first aspect the invention comprises an image processing 

5 . apparatus for the reconstruction of time-dependent representations I(x,t) of an object 
The apparatus may be realised by a computer together with appropriate programs to 
execute the required image processing. The reconstructed representations may 
especially be two-dimensional or three-dimensional cross-sectional images of the 
object, and the object may particularly be the vessel system of a patient where perfusion 

10 shall be investigated or where the flow of a contrast agent shall be reconstructed. As 
perfusion is a time-varying process, the respective cross-sectional images are not 
stationary but time-dependent, too. The apparatus is adapted to approximate the 
function I(x,t) by a given parametric model function I*(a(x),t) and to estimate the 

« 

N-dimensional parameter vector a(x) with the help of a set of projections p'j of the 

1 5 object (2) generated at times tV To accompUsh this functionality, the apparatus 
comprises the following components: 

An approximation module with memory storing the N-dimen- 
sional parameter vector a(x) of the predetermined parametric model 
function I*(a(x),t) that approximates the function I(x,t), wherein I(x,t) 

20 describes the representation of the object Each value of the function 

I(x,t) may for example describe the attenuation coefiBcient or the image 
signal (e.g. the grey-value) of a particular voxel x at time t If 
three-dimensional cross-sectional images are reconstructed, I(x,t) is a 
function depending on four independent variables. The model function 

25 I*(a(x),t), in contrast, depends on the N-dimensional parameter vector a 

and time t in a known manner, because I*(a(x),t) is constructed or pre- 
selected by the user of the apparatus. All you have to know to describe 
the representation of the object (e.g. a cross-sectional image) with 
I*(a(x),t) is then the spatial dependence of the parameter vector a(x). In 

30 other words the problem of finding a spatially dependent functional l(x,0 

for a continuum of times t is reduced to finding N spatially dq>endent 
functions a(x). 
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An input module for the reception of a set of projections p j of the 
object generated at times The input module may for example coiq}le 
the image processing apparatus to an X-ray device generating said 
projections. 

S - An estimator or estimation module that is adapted to estimate the 

N-dimensional parameter vector a(x) with the help of a set of projections 

• * • 

p j of the object, whereia each projection p j is produced at some time tj, 
and wherein the times tj are typically different from each other. Possible 
ways to realize such an estimation module are described in connection 
10 with preferred embodiments of the invention. 

By approximating the real function I(x,t) by a parametric frmction 
I*(a(x),t), the apparatus reduces the complexity of the reconstruction problem 
considerably. Moreover, this approximation has the advantage to allow an estimation of 
the parameter vector a and therefore of the desired representation of the object by 

s 

1 5 projections p j that may in principle be arbitrarily distributed over time. In particular, it 
is not necessary to produce whole sectional images of a (quasi) stationary situation for 
various time points which is for example not possible in perfusion imaging with a 
C-arm based X-ray system. With the described apparatus, however, such systems 
become applicable, and the projections generated by a C-arm system during one or 
20 more sweeps may be used. 

Further, the invention provides an X-ray system suitable for determining a 3D 
dynamic process in an object, the system comprising 

an x-ray source and an x-ray detector placed at opposite positions with respect 
to an examination space and simultaneously rotatable around said examination space 
25 for geuCTating a plurality of x-ray projections; 

a data processing unit for deriving from said plurality of x-ray projections a map 
of the time dependent 3D dynamic process in the object; 

whereby the 3D dynamic process is approximated by a predetermined model 
with a limited set of parameters; 
30 whereby the data processing unit is arranged to estimate parameters in said 

limited set of parameters out of data in the x-ray projections. 
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According to a preferred embodiment of the invention, the apparatus 
comprises an evaluation module for the determination of a perfusion map from the 
representation I(x,t) of a vessel system. If the representation I(x,t) and its approximation 
I*(a(x),t) for instance describe for each time t a cross-sectional image of the object, 
5 perfusion parameters may readily be calculated therefrom and represented in a 
perfusion map. 

The representation I(x,t) of the object and its approximation I*(a(x),t) 
may in general relate to any spatially distributed entity characterizing the object and 
derivable from projections. In a special embodiment that will be considered in more 

10 detail below, the functions I(x,t) and I*(a(x),t) describe for each point t in time a cross- 
sectional image through the object. 

According to a first approach to calculate the desired parameter vector 
a(x) of a function I(a(x),t) representing a cross-sectional image, the apparatus is adapted 
to use an update function AI(x, p*^\ l\x)) for the estimation of the parameter vector a, 

1 5 wherein the update function AI(x, p*^^ is taken from an iterative algorithm for the 
reconstruction of stationary cross-sectional images I(x). One important example for 
such an algorithm is knowii as "Algebraic Reconstruction Technique" (ART). Update 
functions based on this algorithm or derivatives thereof are especially suited for the 
present estimation method. In the i^nlate function, p'^^ is a projection of an object used 

20 in the k-th iteration step, and I^'Cx) is the k-th estimate for I(x). The application of 

update functions of known algorithms has the advantage that well known mathematical 
procedures and existing tools may be exploited. 

In a ftirther development of the apparatus described above, the parameter 
vector a(x) is iteratively approximated by a sequence a^(x), wherein the (k+l)-th 

25 iteration comprises the following steps: 

a) The computation of k-th order estimates I*(aXx),t*j) for at least N 
of the times t j, wherein the times t j may in principle be selected 
arbitrarily- This selection of times may be described by corresponding 
index sets A and B according to ie A and j eB (wherein A and B 

30 normally change from iteration step to iteration step). 

b) The computation of corresponding updates 

AI^'^j = AI(x, p^j, I*(a\x),t*j)), 
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wherein AI is the update function of a known iterative algorithm like 
ART. The update function makes use of said estimates I*(a\x),t j) and the 
measured projections p j that correspond to the selected times t j with ieA, jeB. 

c) The calculation are the new estimate a*^^(x) for the parameter 
5 vector by minimising the following function with respect to 

ieAjeB 

The procedure described above will be discussed in more detail in the 
description of the accompanying figures. Its advantage is that it may be carried out 
iteratively based on standard algorithms like ART and thus may readily be 
10 implemented. 

The projections p'j of the object that are used for the estimation may in 
principle be arbitrarily distributed over projection angle and time, wherein the 
distribution should of course be such that no important developments in the processes 
to be observed in the object, e.g. a perfusion process, are missed. Preferably, however, 

15 the set of measured projections is structured into a number M of subsets, wherein each 
subset comprises only projections p j, j 1> • • •» Q» which were taken from the same or 
approximately the same direction at different times. Moreover, the niunber Q of these 
projections in each subset is preferably greater or equal to the dimensionality of the 
parameter vector a, i.e. Q > N. In other words there are M different projection 

20 directions, and for each of these projection directions a number Q ^ N of single 
projections is produced at different time points. Such a structure of the projections 
allows to use at least N projections taken from the same direction in an iterative 
procedure for the estimation of a(x). Thus, for example, the iterative algorithm 
described above may be carried out such that |A| = 1 in each iteration step, provided 

25 that the times t j with the same superscript i belong to the same projection direction. The 
use of projections from the same direction in each iteration step makes the 
approximation procedure more stable. It should be noted at this point that the separate 
indices i, j of the projections p j were introduced to reflect the aforementioned structure 
of the set of projections; if there is no such structure, the indices i and j are simply 

30 symmetrical, i.e. they play the same role. 

A second approach to calculate the desired parameter vector a(x) for the 
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parametric model function I*(a(x),t) may be based on the optimisation of an objective 
function that compares measured and calculated projections. Especially, the objective 
function may measure the total deviation between all measured projections p j and the 
corresponding projections ft I*(a(x),t j) that are calculated from the model function I*. 
5 Such an objective fimction may for instance be defined by the formula 

ij 

An advantage of this approach is that it may readily be added to the usual 
image reconstruction procedures without a need to change them. Moreover, there are 
a priori no constraints on the kind of the used projections p j; they may for example all 
10 belong to different projection directions and time points. 

According to a further development of the apparatus, the estimation of 
the parameter vector a makes use of an anatomical reference data set This reference 
data set may for example be a stationary cross-sectional image of the object that was 
produced in advance. The use of such a reference data set allows it to reduce the 
1 5 number of independent components in the parameter vector by one, which simplifies 
the estimation process accordingly. Moreover, such a reference data set is especially 
usefiil in combination with the second approach described above where it may provide 
an appropriate starting value for the optimization process. 

The invention further relates to an X-ray examination system comprising 

20 the following components: 

A rotational X-ray apparatus, for example a C-arm based X-ray 
system or a multi-slice CT system, for generating X-ray projections p j of 
an object from different directions. 

An image processing apparatus coupled to the X-ray apparatus 

25 and adapted to approximate the function I(x,t) that exactly describes the 

representation of an object by a given parametric model function 
I*(a(x),t), and to estimate the N-dimensional parameter vector a(x) with 
the help of a set of projections p j of the object, wherein each projection 
p j is produced at some corresponding time t j. The apparatus may 

30 especially be realized by a computer together with appropriate programs 

to execute the required image processing. Moreover, the image 
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processing apparatus may particularly be of the kind described above. 

Because the image processing apparatus allows to reconstruct time- 
dependent representations of an object from projections taken at different stages of the 
processes to be observed in the object, the rotational X-ray apparatus may especially be 
5 a conventional X-ray system with a C-arm. X-ray projections generated during one or 
more sweeps of such a system may be used for the reconstruction of the desired 
intensity function I*(a(x),t). 

According to a further development of the invention, the X-ray 
examination system described above comprises an injection system for injecting 
10 contrast agent into the blood flow of a patient. The injection system may for example 
comprise a catheter and a pumping unit for controUably introducing contrast agent 
through the catheter into the vessel system of a patient The injection system allows to 
perfonn perfusion imaging or reconstruction of contrast agent flow with the 
examination system. 

15 The invention further comprises a method for the reconstruction of time- 

dependent representations of an object comprising the following steps: approximation 
of the function I(x,t) describing the representations by a given parametric model 
function I*(a(x),t); and estimation of the N-dimensional parameter vector a(x) with the 
help of a set of projections p'j of the object produced at times t j. The projections p j may 

20 particularly be generated with a G-ann system or a multi-slice CT system. The method 
comprises in general form the steps that can be executed with an image processing 
apparatus of the kind described above. Therefore, reference is made to the preceding 
description for more information on the details, advantages and improvements of that 
method. 

25 Moreover, the invention comprises a computer program for enabling 

carrying out a method of the kind described above. Such a computer program together 
with a general-purpose computer may especially constitute an image processing 
apparatus of the kind described above. 

Finally, the invention comprises a record carrier on which a computer 

30 program of the kind described above is stored. The record carrier may be any kind of 
permanent or volatile storage medium like RAM, ROM, tord disk, floppy disk, CD, 
magnetic tape or the like. 
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These and other aspects of the invention will be apparent from and 
elucidated with reference to the embodiment(s) described hereinafter.' 



5 In the following the invention is described by, way of example with the 

help of tihie accompanying drawings, in which: 

Figure 1 is a principle sketch of an X-ray examination system and the 
method to generate projections for perfusion imaging according to the present 
invention. 

1 0 Figure 2 is a diagram of a typical course of time-dependent grey-values 

and its approximation according to the present invention. 

The drawing in the \xpper left comer of Figure 1 schematically depicts 
the setup of an X-ray examination system according to the present invention. The 
examination system comprises a rotational X-ray apparatus 3 which is adapted to 
1 5 . generate X-ray projections of an object from different directions d', for example of the 
vessel system 2 of patient to be examined. The rotational X-ray apparatus 3 may 
especially be a conventional C-arm based X-ray system comprising an X-ray source 4 
and an X-ray detector 6 connected via a C-ann 5. The X-ray apparatus 3 is connected to 
an image processing apparatus 1 , which may for example be a workstation with 
20 appropriate software running on it. 

During perfusion imaging or reconstruction of flow in medical 
diagnostics, a contrast agent is injected into the vessel system 2 of patient by a catheter 
(not shown), and the spreading of this injection is monitored in order to gain 
infonnation on the perfusion properties of the vessel system. Due to tiie temporal 
25 development of die perfusion process, die observed image signal of a voxel at 
position X is not constant but usually varies over time. Figure 2 shows a typical 
time-dependence of the grey-value I(x,t) of a voxel x during a perfiision study (black 
line). The complete temporal development in the volume to be observed must tiierefore 
be described by an intensity function I(x,t) which depends on the three-dimensional 
30 position vector x (in cases of tiiree-dimensional image leconstruction) and the 

one-dimensional parameter time t. 

An important step of the procedure to be described in the following is die 
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approximation of said function I(x,t) by a model ftmction I*fe(x),t) that is 
pxedetennined by the user of the system and that depends on an N-dimensional 
parameter vector a and on time t in a known way. The only unknown relation will then 
be the spatial dependence of the parameter vector a(x). The grey curve in Figure 2 
5 shows a typical example of a model function I*(a(x),t) with an 8-dimensional parameter 
vector a(x), wherein pairs of components (aj, ai+i) describe one point in the temporal 
course of the model function and wherein these points are connected by straight lines. 
Of course a lot of other different approximations are suitable, too, for example 
polynomial functions or the like. The user will typically choose a model function such 

10 that it will describe the a priori expected course of intensity values with a minimal 
number of independent parameters. 

In the foUowing, a preferred procedure for the estimation of the 
parameter vector a(x) in the setup described above will be described. The procedure 
comprises the generation of projections p j with the X-ray apparatus 3 from different. 

15 projection directions d' while a perfusion process is going on in the vessel system 2. In 
a preferred embodiment, the X-ray apparatus 3 is rotated forward and backward in such 
a way that for each of M projection directions d' at least N projections p j (1 :^ j < N) at 
different time points t j are acquired. 

From the state of the art iterative reconstruction algorithms such as ART 

20 (Algebraic Reconstruction Technique) are well known (cf. e.g. G. T. Herman, B. 
Meyer; Algebraic Reconstruction Techniques can be made Computationally Efficient; 
IEEE TMI 12 (3), p. 600£F, 1993; F. J. Beekman, C. Kamphuis; Fast Ordered Subset 
Reconstruction for X-ray CT; IEEE Nuclear Science Symposium Conference Record; 
Volume 2, 2000). In these algorithms, the iterative reconstruction of a stationary 3D 

25 image I(x) is achieved with projections p* of an object and the recursion formula 

l'^\x) = l'ix) + AI(x;p^^'\l\xy)^ (1) 

wherein l\x) approximates the unknown function I(x). In many 
algorithms (including ART), AI depends only on the projection of I^'Cx) taken from the 

same direction as p*^\ 

When extending such an algorithm to the non-stationary situation of 
perfusion imaging as it was described above, a similar formula is needed for the 
iterative approximation of the parameter vector a(x) by a sequence a\x). 
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Assume that the k-th element a^(x) of this sequence is akeady known and 
define I^(x,t) := I*(aV)> t). The next iteration step for finding ^\x) must be carried 
out for one of the time points for which a projection p j is available. For this fixed 
time point t*j, the recursion formula (1) may analogously be appUed, which yields 

/**^ (X, r)) = /* (X,/;) + A/(x; (x,^;)) 

5 =r(fl*(x),r;)+A/(x;/7;,r(a*(x),/;) 

From this equation the difference between the model function 
I*(a^^(x),t^j) at the new (imknown) parameter vector and the (known) updated 

function I^^ (x,t j) can be calculated as: 

10 . A reasonable postulation is now to make \xj\ as small as possible by an 

appropriate choice for a*^^(x). As the vector a has N components, this postulation alone 
would however have many solutions and thus convergence cannot be guaranteed. 
Therefore, it is better to adapt several (preferably all) components of a*^\x) 
simultaneously by requiring the minimisation of 

wherein A, B are index sets that may in principle be chosen arbitrarily 
for each iteration step, and wherein the number (|AHBD of the tenns is preferably 

greater of equal to N. 

In order to guarantee a stable convergence (due to a balanced 
20 contribution of all components of a to the minimization), it is preferred to choose for 
the iteration step k -> (k+1) only projections taken from the same direction d\ This 
means that jA] = 1 and |B| ^ N in equation (2). Especially the function 

;ir^(^)= tfcy =i:(/*(5*"(^).^0)-/*(^'W.^0)- A/^'OWJ . 

may be minimized in order to find for each voxel x the new parameter 
25 vector a^^^(x). This solution may be found, for instance, by a least-squares method or 
some other suitable procedure known from numerical mathematics. 
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In sununary, each iteration step of the reconstruction algorithm can be 

realized by 

Computation of N estimates I*(a^(x),t^j) for the time points 
^i,„.,t*N based on the already known parameter vector a^(x) and an 
5 (arbitrarily chosen) index i. 

Computation of N updates A/*-*!,..., A/^-'jv using the acquired 

« 

projections p j and said N estimates. 

Computing during backprojection for each voxel x the updated 

parameters a^'^^Cx) by minimizing z^M- 
10 As variants, the algorithm can also be carried out if the projection 

directions of the N projections pS,.,.,p^N differ slightly. In particular, the projections 
acquired with an angular difference about 1 80° at the turning point of the C-arm 
movement can be considered simultaneously within a reconstruction iteration to 
improve stability. Furthermore, an anatomical reference data set can be used to reduce 
15 the number N of parameters (jc),..., a^v (x) by 1 . 

According to another approach, the parameter vector a(x) may be found 
by minimisation of an objective function like 

wherein 1 ^ i ^ M, 1 ^ j ^ N, and wherein Pi denotes the operator 

■ 

20 projecting the image I*(a(x),t) according to the projection geometry d'. In this case it is 
not necessary (or helpful) that all projections p'j with the same superscript i belong to 
the same projection direction. In contrast to the diagram of Figure 1, the p j may 
therefore be quite arbitrarily distributed over projection angles and time. 

Mmimisation of the function %^ may e.g. be achieved by an iterative 

25 procedure. In this case it is preferred to use a (2D or 3D) sectional image of the region 
of interest as a starting value for I*(a(x),t*j). Such a "starting image" may be produced 
by any appropriate method, e.g. by backprojection from a set of projections. Moreover, 
the starting image may be taken from previously (e.g.. days before the perfusion 
imaging) recorded data, thus minimising the X-ray exposure of the patient Preferably, 

30 however, the starting image is produced immediately before the projections p j in order 
to correspond to the situation on these projections as close as possible. 
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The invention enables perfusion imaging on C-arm based X-ray systems 
as it does not require acquisition of (multiple) cross-sectional images with a frame rate 
in the order of 1 second. When applied on multi-slice CT (MSCT) systems, the 
invention facilitates perfusion protocols with considerably reduced dose. 
5 Finally it is pointed out that in the present application the term "compris- 

ing" does not exclude other elements or steps, that "a" or "an" does not exclude a plura- 
lity, and that a single processor or other unit may fulfill the functions of several means. 
Moreover, reference signs in the claims shall not be construed as limiting their scope. 



